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^ ! Abstract 
> 

in 

^ I The vacuum polarization due to a virtual heavy quark pair and specif- 

^ ' ically the coefficients of its Taylor expansion in the external momen- 

^ ' turn are closely related to moments of the cross section for quark- 

CD ' antiquark pair production in electron-positron annihilation. Relating 

measurement and theoretically calculated Taylor coefficients, an ac- 
curate value for charm- and bottom-quark mass can be derived, once 
corrections from perturbative QCD are sufficiently well under con- 
■ trol. Up to three-loop order these have been evaluated previously. We 

now present a subset of four-loop contributions to the lowest two mo- 
^ ■ ments, namely those from diagrams which involve two internal loops 

• from massive and massless fermions coupled to virtual gluons, hence 

of order a^rij. The calculation demonstrates the applicability of La- 
porta's algorithm to four-loop vacuum diagrams with both massive 
and massless propagators and should be considered a ffist step to- 
wards the full evaluation of the order contribution. 
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117312, Russia. 



1 Introduction 



The correlator of two currents is central for many theoretical and phenomeno- 
logical investigations in Quantum Field Theory (QFT) (for a detailed review 
see e.g. Important physical observables like the cross section of electron- 

positron annihilation into hadrons and the decay rate of the Z-boson are related 
to the vector and axial-vector current correlators. Total decay rates of CP even 
or CP odd Higgs bosons can be obtained by considering the scalar and pseudo- 
scalar current densities, respectively. Some of these studies require to calculate 
the correlator for arbitrary momentum q^. For many applications, however, the 
knowledge of a few derivatives at = is sufficient. 

Two-point correlators have been studied in great detail in the framework of 
perturbative QFT. Indeed, due to the simple kinematics (only one external mo- 
mentum) even multi-loop calculations can be performed analytically. The results 
for all physically interesting diagonal and non-diagonal correlators (vector, axial- 
vector, scalar and pseudo-scalar) are available up to order a^, taking into account 
the full quark mass and momentum dependence [SI IH E] ■ 

The determination of the heavy quark masses with the help of QCD sum rules 
requires the detailed knowledge of the heavy quark correlators. In fact, as noted in 
ref . jS] , the precise determination of the charm- and bottom-quark masses would 
be further improved by including four-loop corrections, hence (9(q;^), to at least 
the few lowest Taylor coefficients of the polarization function. 

Technically speaking, the moments i.e. certain weighted integrals of the cross 
section for electron-positron annihilation into heavy quarks, can be expressed 
through massive tadpoles or vacuum diagrams (diagrams without dependence on 
the external momentum). The evaluation of these "massive tadpoles" in three- 
loop approximation has been pioneered in [7j and automated and applied to a 
large class of problems in [8, . However, in spite of huge progress in calculational 
techniques during recent years the problem of analytical calculation of massive 
tadpoles at the four-loop level has not yet been mastered.^ 

Similar to the three-loop case, the analytical evaluation of four-loop tadpole 
integrals is based on the Integration- by-parts (IBP) method ||9j. In contrast to 
the three-loop case the manual construction of algorithms to reduce arbitrary 
diagrams to a few master integrals is replaced by a mechanical solution of a host 
of again mechanically generated IBP equations fTH] . 

Unfortunately, the price for this automatization is an enormous demand on 
computational power. A system of more than ten million linear equations has 

^The numerical evaluation is certainly possible for individual contributions; however, one 
could hardly imagine a direct numerical evaluation of hundreds of thousands of separate terms 
which appear after performing necessary expansions and traces at the four-loop level. 
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to be generated and solved. In the present publication we present only a partial 
result. We restrict ourselves to the first two non- vanishing moments and consider 
only four-loop diagrams with the maximal number (two) of closed fermion loops 
inside. This corresponds to terms proportional to nf, n| and niUh, where n; 
denotes the number of light quarks, considered as massless and the number of 
massive ones. This leads to a system of about one million equations. 

In general, the tadpole diagrams encountered during our calculation contain 
both massive and massless lines. As is well-known, the computation of the four- 
loop /9-functions requires the consideration of four-loop tadpoles only composed 
of completely massive propagators. Calculations for this particular case have been 
performed in ^TJ ^] . 

The outline of this paper is as follows. In section |2l we briefly introduce the 
notation and discuss generalities. In section IHl we discuss the reduction to master 
integrals, describe the solution of the hnear system of equations and give the 
result for the 0{alnj) contribution to the polarization function for the lowest 
two moments. Our conclusions and a brief summary are given in section |3] 



2 Notation and Generalities 

The vacuum polarization tensor 11^*^ (g^) is defined as 

W\q') = z J dxe'^^0\Tf{x)f{0)\0), (1) 

where q'^ is the external momentum and is the electromagnetic current. The 
tensor n^'^(g^) can be expressed by a scalar function, the vacuum polarization 
function n(g^) through 

U^-{q^) = ^ n(g2) + q^ n^(g2). (2) 

The longitudinal part Uiiq^) is equal zero due to the Ward identity. 
The confirmation of 11^^ (g^ = 0) = and 11/, (g^) = will constitute an important 
check of our calculation. The constant n(g^ = 0) relates the QED coupling in the 
on-shell scheme and the MS-scheme. The first and higher derivatives of n(g^) at 
q^ = contain important scheme independent information and will sometimes be 
called the physical moments. The imaginary part of the polarization function is 
related to the physical observable R{s) 

R{s) = 12 7rlmU{q'^ = s + ie). (3) 

and properly weighted integrals of R{s) obviously coincide with the Taylor coef- 
ficients 
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This justifies to study the Taylor expansion around = 0. In this case the 
coefficients are given by massive tadpole integrals 

= ^.i:CnZ- (5) 

n>0 

with the dimensionless variable z = g^/(4m^), where m is the mass of the heavy 
quark and Nc denotes the number of colors. It is convenient to define the ex- 
pansion of the coefficients C„ of the polarization function in the strong coupling 
constant as 

a = Ci'^ + a. + al C^^ + alC^^ + ... (6) 

with a, = — . Within this work we consider the contribution of Cp-* and 
define C^'>\^2 = T"^ Cp C^f\ where T denotes the normalization factor of the 
fundamental-representation generators defined by Tr[)f:'^t''] = T 5""^ and is 
the Casimir operator in the fundamental representation. 



3 Calculations and Results 
Reduction to master integrals 

In the ffist step the reducible scalar products in the numerator of the integrands 
have been removed in the sense that trivial tensor reduction has been performed. 
All reducible scalar products have been expressed in terms of their associated 
denominators. Through this the remaining integrals can be mapped upon a set 
of 12 independent topologies. 

Through the expansion in the external momentum q the derivatives acting on 
the polarization function generate additional powers of the denominators of the 
integrands. The deeper the expansion is, the higher powers are obtained. Let Md 
denote the sum of powers of propagators Di minus the number of propagators 
of the generic integral. Let Mp denote the total sum of powers of the irreducible 
scalar products in the numerator of the integrand. Then one obtains integrands 
with Md up to 8 and Mp up to 3, for an expansion of the n^-contribution up to 
the first physical moment. 

For the reduction of these integrals to a set of a few master integrals the 
standard method of IBP has been used. The reduction was implemented following 
the ideas described in ref. jTUl CSl El • In order to reduce the polarization function 
to master integrals for an expansion up to the first physical moment a system of 
around 880000 equations has to be generated and solved. For a deeper expansion 
of the polarization function one obtains higher values for Md and Mp. This requires 
a huger system of IBP-identities in order to obtain a reduction to master integrals. 
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Figure 1: List of independent topologies belonging to the nj-contribution. The 
topologies Ri-Rq are reducible, whereas the topologies Mi-Mg are master inte- 
grals. The wavy and solid lines depict massless and massive propagators respec- 
tively. 

A lexicographical ordering has been introduced assigning to each integral a weight 
describing its "complexity" . Integrals with increasing powers of the denominator 
and increasing number of irreducible scalar products are denoted as increasingly 
complicated. The linear system of equations has been solved with a program based 
on FORMS jT^I which uses FERMAT |TB] for simplifying the rational functions in 
the space-time dimension which arise in this procedure. Complicated integrals 
are systematically expressed in terms of simpler ones and then substituted into the 
other equations. Some contributions have been checked independently with the 
program SOLVE J7j from which also some experience concerning the procedure of 
ordering the integrals according to their complexity has been gained. Masking of 
large integral coefficients is used, a strategy also adopted in the program AIR [TH] . 

Exploiting the symmetries of the diagrams by reshuffling the powers of the 
propagators of a given topology in a unique way strongly reduces the size of the 
initial input and, similarly, in the second step the number of equations which 
need to be solved. 

The solution of the system leads to a set of a around 130000 independent 
equations. With the help of these equations the first two Taylor coefficients of the 
polarization function can be expressed in terms of 6 master integrals Mi-Mg, with 
denominator powers one and no irreducible scalar products (M^ = 0, Mp = 0), 
which are shown in Fig. H 
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Calculation of the master integrals 

The master integrals which belong to the diagrams Mi — Mg are defined 
4 — 2e space-time dimensions through 



in d 
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and the normalization factor 



d% 1 



{2t:YDI 16 7r2 



V{l + E) 



m 



47r/i2 



(13) 



(14) 



The factor fi denotes the renormalization scale. 

Before calculating the master integral Mi we consider at first the following 
combination of integrals with dots, where a dot on a line denotes an additional 
power of the associated denominator 



m 




iVio + 0{e) 



(15) 



The three-loop topology in the left hand side of eq. ()15|) as well as all fol- 
lowing three- loop diagrams are normalized by (of eq. ()14|) ). The combi- 
nation in eq. ()15|) is finite and can be integrated numerically, with the result 
A^io = 5.3111546.... 
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The relation between the dotted topologies in eq. (fT^ and master integrals can 
be obtained via IBP. One finds the following relations: 





d-sy* (2d- 7) 
16 {d-4) 



3(d-3) (3 d-10) (3d-8) 
256 (d-4) m6 



+ 




-3) (d-2) (3d-8) (7d^-48 d+82) / 
128(d-4) (2d-7)m8 



{d-2f (lld-38) 
128 (d-4) m» 





(16) 



(17) 



The three-loop integrals in the right hand side of eq. (fT7j) and the factorizable 
four-loop amplitude in the right hand side of eq. ()l(i|l can be calculated with 
MATAD |i8j or taken from ref. Inserting eq. ()16p and (fT7|) into eq. (fT^ 
leads to^ 
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- — (4335 - 1440 a4 - 60 log*(2) + 120 A^io + 60 log^(2) 
15 



TT 



with 
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00 00 
Cn = E ^ and 04 = E 



2^fc^ 



:i8) 



(19) 



fc=i fc=i 
The same procedure has been applied for calculating the master integral M2. The 
topology M2 with 3 symmetrical distributed dots is finite 



m 




N20 + 0{e) 



(20) 



and can be integrated numerically. The calculation yields A''2o = 5.40925606 . 
The master integral M2 can be calculated using the IBP identity 




+ 



((j-4) (d-3)^ (2d- 7) 
6(3d-13) (3 d-ll)m6 



(d-2)^ (d-1) 
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(21) 



inserting eq. (PT|l into eq. (pUjl and solving with respect to M2 results in: 
3 



Mo 



m 



2e' 



19 67 1 \ 

+ ^ + ^ + ^(127-6iV20 + 21C3) + O(^°)). 



^We have been informed that the same result has been independently obtained in PU] , 



6 



This result is in agreement with eq. (4) in ref. 
The calculation of the master integral M3 is easy, 

^3 = ^'fA + A + A-7^(39- 38C3)-7^(4710 + 23 7^^ 
\6€^ Ge'^ 2e^ 6e 60 

- 1900 C3) - — (5934 + 23 TT^ - 684 Cs - 2580 Cs) + 0{e^)] . (22) 
For completeness we also give the results for the factorized master integrals 

- 184320 04 - 7680 log^(2) + 7680 log2(2) tt^ + 2176 vr^ - 228480 Cs) + 0{e^) 

/ 1 3 43 1 1 
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V 3e^ 2e3 1252 24e ^ ^'^^ 240 ^ 



2880 Cs) + (60435 + 288 vr^ - 13760 Cs - 23040 C5) + 0{£^' 
480 ^ 



X , 1 4 10 20 
Me = m« (- + - + — + — + 35 + 56£ + 0(e")) . 



Result for the 0{cx^'ny) contribution 

Inserting the above master integrals into the reduced nj-contribution and per- 
forming renormalization of the strong coupling constant a^, the external current 
and the mass m = m(/i^) in the MS-scheme, leads to the following result for 
the first two moments of the four-loop (^(afnj) contribution of the heavy quark 
vacuum polarization function: 

.(3) _ . . ^7043 127 ^ ,1^/ 37 7 N , 2 . 2 ,3^ 
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+ 1^23328 54^' 81 + 27 > 27 ^ 
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/ 3779 203 X 472 16 
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2 /42173 _ 112 . _ 1784 236 ,2 _ 8 A 
^ V32805 ~T35^'~ 3645 ^^^ 3645 '^"TSS^^^ ^ ^ 

with = log(^). 

The n^^-contribution has been checked independently by taking into account the 
corresponding two-loop case, in which the gluon propagator has been replaced 
by a gluon propagator containing a renormalon chain with two massless fermion 
one-loop insertions. The computation has been performed in a general ^-gauge 
and it has been checked that the dependence on the gauge parameter vanishes. 
Furthermore it has been checked that both coefficient functions Cq and Ci meet 
the standard renormalization group equation. Numerically one finds for the coef- 
ficient Co and Ci at /i — m: 

Co = a, 1.4444 + (1.5863 + 0.1387n/j + 0.3714n;) (25) 

+ al (0.0252 Uh m + 0.0257 - 0.0309 , 

Ci = 1.0667 + as 2.5547 + (0.2461 + 0.2637 Uh + 0.6623 ni) (26) 

+ al (0.1658 Hh ni + 0.0961 + 0.0130 nl) . 



4 Summary and Conclusion 

Using the IBP method and Laporta's algorithm, we have evaluated a gauge in- 
variant subset of the four-loop tadpole amplitudes contributing to derivatives of 
the vacuum polarization at — 0. All loop integrals have been mapped on a 
minimal set of independent topologies. Then an elaborate automated procedure 
has been developed and applied, which identifies equivalent amplitudes, factor- 
izable contributions, discards massless tadpoles, and performs symmetrization. 
Solving a system of nearly one million linear equations, all amplitudes can be 
expressed through six master integrals. These have been evaluated analytically 
or numerically to high precision. The present work can be seen as a first step 
towards the evaluation of the full set of four-loop amplitude contributions to the 
vacuum polarization function. 
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